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ABSTRACT OF THE DISSERTATION 



Comparison of Momentum, 
Sensible and Latent: Heat Fluxes 
Over the Open Ocean 
Determined by the Direct Covariance, 
Inertial and Direct Dissipation Technique 



by 



Gregory Frank Dreyer 

s 

Doctor of Philosophy in Oceanography 
• University of California, San Diego, 1974 
Professor Carl H. Gibson, Chairman 

Measurements of velocity, temperature and humidity were made 
in the atmospheric surface layer over the open ocean from the 
Research Platform FLIP, to compare estimates of the fluxes of 
momentum, sensible and latent heat, using the direct covariance, 
inertial and direct dissipation techniques. Measurements were made 
during a series of oceanographic cruises approximately 60 miles off 
the coast of Baja California, between August 1970 and April 1973. 
Direct covariance estimates of the fluxes were corrected for 
instrument tilt using simultaneous measurements of pitch angle, and 
dissipation estimates were corrected for stability effects. 

Agreement was obtained between estimates of momentum 



XXIV 



flux by the dissipation techniques and the direct measurements within 
the uncertainties of the direct estimates (± 25%) , applying the 
assumption that total production of kinetic energy (mechanical + 
buoyant) equals dissipation. From the direct dissipation technique, the 
average value determined for the velocity subrange constant was 

= 0. 53 ± 0. 03 , Trends in the direct dissipation estimates indicate 
that vertical turbulent transport and pressure transport of kinetic 
energy cannot be neglected and may be approximately half of the 
estimates of dimensionless transport measured over land. Comparisons 
were made between directly measured momentum fluxes corrected for 
instrument tilt using measured pitch angles, and values obtained using 
estimates of pitch angle in lieu of measurement. Pvesults of the 
comparison suggest that large errors in direct momentum flux 
estimates can be incurred if tilt corrections are based on estimates of 
pitch angle. 

For agreement between inertial dissipation estimates and 
directly measured latent heat fluxes, an average value of the humidity 
subrange constant = 0.21 i 0. 05 was required. The direct 

| q 

dissipation technique was not applied to estimate latent heat fluxes since 
the frequency response of the Lyman-alpha humidiometer was not 
sufficient for direct estimates of the dissipation of humidity variance. 
Dimensional humidity spectra were similar to horizontal velocity 
spectra, and water vapor flux cospectra were similar to momentum 
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^flux cospectra. 

To- obtain agreement between direct estimates of sensible heat 

flux and inertial dissipation estimates, an average value of the 

temperature subrange constant B = 2. 2 ± 0. 3 was required. 

0 

Estimates of sensible heat flux by the direct dissipation technique were 

higher than direct measurements by a factor of 2 or more. From 

direct estimates of the dissipation of temperature variance (± 20%) 

and sensible heat flux (± 10%) , dissipation was greater than production 

of temperature variance by as much as a factor of 4 , not accounted 

for by uncertainties in the direct estimates, or by temperature flux 

s 

divergence. Direct dissipation estimates of /? varied from 0. 7 i 

0 

0.07 to 2. 4 ± 0. 2 . Temperature spectra from different cruises with 
different environmental conditions did not compare with one another, 
and were not similar to humidity or velocity spectra. Cospectra of 
sensible heat flux exhibited similar differences. The differences 
between production and dissipation of temperature variance, and the 
behavior of the temperature spectra and sensible heat flux cospectra 
could be accounted for by consideration of sources of production usually 
neglected in the temperature variance budget, and attributed to the 
combined effects of ocean spray evaporation and radiative heating/ 
cooling. 

A wide range of bulk coefficients were obtained from comparison 
of estimates of sensible heat flux by the direct covariance and bulk 
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aerodynamic techniques, indicating the bulk technique is not reliable 
for estimates of sensible heat fluxes. A value of (1. 29 ± 0. 36) x 10 
was obtained for the latent heat flux bulk coefficient. 
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1. INTRODUCTION 



Energy and moisture transfer from the oceans to the marine 

atmosphere are major driving faetors of atmospheric and oceanic 

eireulation. Within the last deeade a vast amount of scientific 

investigation and experimentation has been directed toward an 

understanding of the energy exchanges and physieal processes which 

oceur at the atmosphere - oeean interface. Aeeurate determination of 

the vertical fluxes of momentum, sensible and latent heat is of prime 

importance. This researeh has as its primary objectives the direet 

\ 

measurement of the turbulent fluxes and comparison of various 
techniques of estimating the fluxes indireetly from related statistical 
quantities. 

Perhaps one of the earliest attempts to estimate the transfer of 
momentum and heat from the oeean was made by G. I. Taylor in 1913 

from a whaling ship using balloon and kite observations of mean 

1 

temperatures and wind velocities (Taylor > 1970). Due to severe 

instrument limitations, early methods of estimating the turbulent 

fluxes have consisted of semi-empirieal theories relating mean 

properties or gradients to the aetual fluxes. No direct flux calculation 

was possible. Many of these methods including the bulk aerodynamie 

method, integral method, and profile method are diseussed by Roll 

2 

(1965). 



1 



With the evolution and development of sophisticated instrumen- 
tation and experimental platforms for use on the open ocean, accurate 
measurements of atmospheric and oceanographic parameters of scales 
ranging from millimeters to kilometers have become possible. As an 

outgrowth of the recommendations of the Joint Panel on Sea-Air 

3, 4 

Interaction of the National Academy of Sciences (1962, 1966) the 

Bardados Oceanographic Meteorological Experiment (BOMEX) was 

conducted during the summer of 1969. The principal objective in air- 

sea interaction was the study of the fluxes of momentum, sensible and 

< 5 

latent heat described by Davfdson (1968) and Kuettner and Holland 

£ 

(1969). Current results from" BOMEX have been summarized by 

7 

Holland (1972). Results of direct flux measurements from the 

research vessel FLIP obtained using the covariance or eddy correlation 

method and the well recognized difficulties of the measurement are 

8 9 

discussed by Pond et al (1971), Phelps (1971), Phelps and Pond 



12 



(1971)^0 and Leavitt (1973). ^ The dissipation technique for 

estimating momentum flux was used by Gibson and Williams (1969) 

from FLIP in 1968, and for sensible heat flux during BOMEX by 

13 

Gibson, Stegen and Williams (1970) and Stegen, Gibson and Friehe 
14 

(1973). 



Experimental data for this thesis was obtained from measure- 
ments made from the Research Platform FLIP during a series of 
oceanographic cruises between August 1970 and April 1973 off the 



coast of Southern California. Mean and fluctuating component r 

horizontal and vertical velocity, temperature, humidity, instil nt 

motion, and mean sea surface temperature were measured 

simultaneously. Fluxes of momentum, sensible and latent heat were 

obtained by the direct eddy correlation technjq\ie corrected for 

instrument /platform motions and compared with estimates of fluxes by 

the direct and inertial dissipation techniques. The influence of 

atmospheric stability on the determination of the fluxes by the latter 

two methods was determined. Latent and sensible heat fluxes were 

also determined using the bulk aerodynamic technique (including 

\ 

stability effects) and compared to estimates of the turbulent heat fluxes 
by the direct covariance technique. 



Z. FLUX TECHNIQUES 



2 . 1 Direct Covariance Technique 

The optimum method of determining the turbulent fluxes of 
momentum and sensible and latent heat is to measure directly and 
appropriately average the covariances of vertical velocity with the 
suitable variable of interest. This method is known as the direct 
covariance or eddy correlation technique. The technique. is based 
on the definition of the vertical fluxes as given by 

Momentum Flux 

\ 

Sensible Heat Flux 
% Latent Heat Flux 

where u , w are the horizontal and vertical fluctuating components 

of the streamwise velocity vector U , 0 is the fluctuating component 

of temperature (° C) , and q the absolute humidity fluctuation 
3 

(jigm/cm ). The constants are p density, C specific heat 

capacity, and L latent heat of vaporization. Overbars indicate 
v 

time averages assuming stationary flow. The covariances are 

obtained by integration of the appropriate cospectrum $ (where x 

is u , 0 , or q) between low and high frequency limits in the same 

manner as the variances of the individual turbulent components are 

obtained from the appropriate spectrum ^ . The frequency range 

xx 

is determined to include all significant contributions to the integrals 



r = - p uw 

H c = pC w6 (1) 

S p 

H„ = L wq 
jC v 



4 



g 



iven by 




$ (f) df and x 2 

wx 




<£ (f) df 

xx 



( 2 ) 



The low frequency (f^) or high pass cut-off is normally determined 

by removal of the record mean on data processing and is thereby set 

by the length of the data record. The low pass cut-off (f ) is 

b 

initially limited by the bandwidth limitations of the instruments used 

to measure the parameters u , w , 0 and q , and is subsequently set 

by low pass filtering during data analysis, usually adequate to 

determine the full value of the appropriate flux* Filters are used to 

reduce high frequency noise and minimize aliasing in analog to digital 

conversion* In this work the frequency range of the integration was 

-3 

between about 1 x 10 and 10 Hz , 

Determination of the fluxes directly by this technique is fairly 
difficult in practice due to the effects of instrument platform motion 

g 

on the measured turbulent velocity components. Pond et al (1971) 
obtained values of the covariance uw with the eddy correlation 
technique using a sonic anemometer to measure horizontal and vertical 

velocity from FLIP during Operation BOMEX. Employing an empirical 

# 

correction procedure, significant corrections to uw were required to 

account for FLIP motion and effects on the flow field. Rotations of 

the principal axis of the Reynolds flux tensor were selected to make 

1 /2 

the spectral correlation coefficient R (f) = $ / ( $ $ ) 

^ uw uw uu v/w 



''Nsqual to a value of -0. 5 for 0. 01 < £ z'/U < 0. 1 as suggested by 
Smith (1967) and Weiler and Burling (]967). ^ Pond et al^ found 
typical correction values for uw were 13% per degree tilt and 5% 
and 3% per degree for wq and w0 respectively, and it was felt that 
the measurements were within 1-2° of the correct coordinate system. 
Contributions to the u , w spectra and uw eospeetra by wave 
induced FLIP motion were not included in the integration when 
determining variances and covariances, Despite the fact that no 
theoretical assumptions or approximations are required to employ the 

g 

direct covariance technique, Pond et al assumed that there was no 
distortion of the Reynolds stress due to the large tilt angles calculated 

10° ) which would affect the correlation between u and w . 

17 

Kaimal and Haugen (1969) have found that the correlation 

coefficient r = ■■ ■ - , varies over a wide range from near zero to 

uw a o 

u w 

-0. 5 , particularly under unstable conditions. This suggests that the 

spectral correlation coefficient R^ w may also vary with stability. 

1 17 

Based on results from the Kansas experiment, Kaimal and Haugen 
suggest a need for ± 0. 1° accuracy in the internal alignment and 
mounting of sensors used to measure the direct covariance (uw) in 
the atmospheric boundary layer, * 

For this research measurements of instrument platform tilt in 
the u , w plane were made using a vertical gyro located at a known 



angle with respect to the velocity sensors. Outputs of the vertical 
gyro were recorded simultaneously with the fluctuating velocity 
components and temperature and humidity fluctuations to determine 
corrected values of the covariances. Description of the instrument 
package and correction procedures are contained in the 
instrumentation and data reduction sections. 



2. 2 Dissipation Techniques 

Both the inertial and direct dissipation techniques are somewhat 
advantageous over the direct covariance method in that they are much 

S, 

less susceptible to errors caused by instrument platform motion. This 
is due to the fact that the flux estimates are based on measurements of 
the small scale structure of the velocity and scalar fields. Three 
variations of the dissipation method differ in the technique used to 
determine the rates of dissipation of velocity or scalar variance are 
determined. The dissipation may be measured directly or inferred 
from the inertial subrange of the appropriate spectra or of the 
structure function. 

The technique for relating the dissipation of velocity and scalar 

variance to the fluxes of momentum and heat was first suggested by 

1 8 

R. J. Taylor (1961). This technique is based on the assumptions 

19 

discussed in Lumley and Panofsky (1964) that mechanical production 
is equal to energy dissipation at the same height in the constant flux 



layer for near neutral conditions, so that 



(production) - uw 



dU 

dz 



€ (dissipation) 



(3) 



where £ = mean rate of viscous dissipation = 2 v e 

ij 

(u , +u, .)/ 2 ; v = kinematic viscosity; p density, 
1 > J h 1 

streamwise velocity, z = vertical height above the m 




U 

ean 



and u and w are defined as previously. 






aiid e s 
k) 

m ean 



sea surface, 



The same assumption is made for the production and dissipation 
of scalar variance, i* e. , 



— V 

— dy y 

•(production) -yw — = 



(dissipation) 



(4) 



where y and y represent the mean and fluctuating scalars (tempera- 
ture (T+0 ) and humidity (q+q)), X “ mean rate of dissipation of 

. ~2 

scalar variance = jgrad y| , and , the molecular diffusivity 

of the scalar. 



Assuming a logarithmic profile for velocity and the scalars, and 
implicitly invoking conditions of stationarity and horizontal 
homogeneity, we find 



dU 

dz 



u^ 

Kz 




= -uw 



(5) 




dz Kz a 

y 



y * 




( 6 ) 



''Nvhere u_, ; is the friction velocity, K is von Karman's constant, 

normally assumed to be 0.4 , y , the scalar scale, and 

Oi v = (yw bU/b z)/(uw by /b z) = K /fv the ratio of eddv transfer 
/ H m 

coefficients for the scalars and momentum, assumed unity for neutral 

conditions (Reynolds analogy). Also assumed is K ^ K for 

H E 

humidity. Expressions for the fluxes can be obtained by combining 
equations (3) and (5), and (4) and (6) 



2 



u, 



(K e 




(7) 




a x 

-. r . y 

2u,K 




( 8 ) 



2. 2a Inertial-dis sipation technique 

J g 

•Taylor determined values of £ and x from the inertial 

0 



subrange of the second order structure functions for velocity and 
temperature. From the similarity hypothesis of Kolmogorov (1941) 



Z0 
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for turbulent velocity fields, and Obukhov (1949) and Corrsin (1951) 
for turbulent scalar fields assuming local isotropy, the second order 
structure functions in the inertial subrange are given by 



2 _ 2/3 2/3 

D (r) = [ U (x+r ) -U(x)] = a € r 

uu • u 

2 -1/3 2/3 

D 0e ( r ) =[T(x+r) -T(x)] = b Q X Q e r 



(9) 

( 10 ) 



where U and T (or humidity) are measured at streamwise separation 



distance r , between L q (energy scale) and n (dissipation scale). 

3-1/4 

7? is also known as the Kolmogorov scale, defined as (v /€ ) 

Values of € and \ can also be determined in the same 
manner from the inertial subrange of the one-dimensional spectra 
given by 



, , - 2/3 -5/3 

$ (k ) = OC C k 

UU U 



( 11 ) 



$ (k) = /g X € “ 1/3 k - 5 / 3 

YY y y 



( 12 ) 



where OL and B are universal constants assumed known, and the 
u y 

\ — 

radian wavenumber k = 2TTf/U by Taylor’s hypothesis. This 

technique requires accuiuite values of the constants a. and (3 (or 

u y 

a and b ) , and the velocity and scalar spectra and structure 
u y 

-5/3 2/3 

functions to follow the k~ (r ) form of the inertial subranges. 

In this work, examples of temperature spectra obtained from 
open ocean atmospheric data are presented which appear to exhibit a 
very limited, if any, -5/3 inertial subrange. Application of Eq.(12) 
to single frequency (wavenumber) measurements of the spectral 
function could lead to large errors in the estimation of sensible heat 
flux if an inertial subrange does not exist. Similar care must be taken 

t 

when applying Eq. (10) to single separation measurements of the 
structure function. A substantial portion of the temperature spectrum 



should be measured (to ~ 100 Hz) to determine if a subrange does exist 



before this method is applied. Humidity spectra appear to exhibit a 
- 5/3 inertial subrange behavior over a wider range of frequencies 
(to ~20 Hz) than temperature 



2* 2b Direct dissipation technique 

This technique involves measurement of the mean square 

values of time derivatives of fluctuating streamwise velocity and 

temperature to calculate dissipation rates directly. Direct estimates 

of £ and X in the atmospheric boundary layer over the ocean from 

0 

FLIP in 1968 and during Project BOMEX are described by Gibson and 

12 13 14 

Williams (1969), Gibson et al (1970) and Stegen et al (1973). ' 

Assuming local isotropy and Taylor’s hypothesis (U d/dx = -d/dt), 

the viscous dissipation may be determined from the relations 



and 



€ 



15 v 








(13) 



(14) 



In practice the appropriate time derivative spectra are integrated to 
obtain values of the mean square time derivatives. Difficulties in this 
technique involve the high spatial resolution required of the sensing 
probes to nearly the Kolmogorov scale ( 77 ) , of the order of 1 mm 
in atmospheric flows. High frequency response (~ 2000 Hz) and very 



high signal-to-noise ratios are required in the sensor and associated 
circuitry. At the present writing a humidity sensor with sufficient 
frequency response was not available to measure adequately mean 
square humidity time derivatives to calculate the humidity dissipation. 
As a result, the direct dissipation technique is applied only to 
estimates of the momentum and sensible heat fluxes. 

The direct dissipation technique does offer several 
advantages over the other techniques. Estimates of the fluxes require 
relatively simple instrumentation and data analysis, and are- not 
affected by sensor motion because of the high frequency response and 
fine scale spatial resolution of the sensors used. Measurement of 
spectra over the entire range of frequencies (wavenumbers) are 
readily obtainable to determine spectral shapes. The spectra 
combined with dissipation measurements provide a means of accurately 
determining the universal constants. 

2. 3 Stability Considerations 

Application of the dissipation techniques using Eqs. (7) and (8) . 

does not account for effects of stability on the flux estimates. These 

equations are modified by consideration of the turbulent budget 

equations for kinetic energy and scalar variance, aiid the flux-profile 

relationships in the surface layer, within the framework of the Monin- 

.23 

Obukhov surface-layer similarity theory (Obukhov, 1946 and Monin- 
24 



Obukhov, 1954). 



In the first tens of meters of the atmospheric boundary layer, 



the fluxes of momentum and heat are assumed constant. Recent 

experimental support for this assumption has been given by Haugen 
Z 5 Z 6 

et al (1971) and Dyer and Hicks (197Z). In this n constant flux 

layer, 11 according to similarity theory, the turbulence structure is 

determined by the surface shear stress T , the surface heat flux H , 

o s 

the buoyancy parameter g/T, and z the vertical height. The latent 
heat flux H is included to account for humidity effects on buoyancy 
in the boundary layer over the ocean. From these pcirameters 
characteristic scales for velocity (u^) , and scalars (}/ = T o , q , ) 
are defined (as in Eqs. (5) and (6)). The length scales are z and the 
Monin-Obukho v length, defined as 



u. 3 T 



L = - 



(15) 



K g w 0 



v 



where the mean and fluctuating virtual temperature, T and 0 , as 

v v 

19 

given by Lumley and Panofsky (1964), include humidity contributions 
to buoyancy. They are written as 



T = T (1. 0 + 0. 61 M) 
v 



( 16 ) 



and 



v 



0 + 0. 61 Tm 



(17) 



where M and m are mean and fluctuating specific humidity (gm/gm). 
The scale L is a key independent variable of similarity theory and 



\ determines the thickness of the surface' layer above which buoyancy 
factors are not significant. 

Applying the scaling parameters to the vertical gradients of 
velocity and the scalars (temperature, humidity), Eqs. (5) and (6) are 
rewritten as 



dU 
d z 



<Lz = 

d z 



u 



77~~ $ (z/L) 

K z m 



V* 

: $ (z/L) 

Kz a' H 

y 



(IB) 



(19) 



\ 

The functions $ , the dimensionless wind shear, and , 

m H 

the dimensionless scalar gradient, are evaluated empirically. It is 

generally assumed that <3? = <£> for temperature and humidity. 

H E 

2 7 

Monin and Yaglom (1971) conclude that the forms of $ and <£ 

m H 

for unstable conditions (the prevalent condition over the open ocean) 
are best given by 



-1/4 

$ = (1 - a z/L) f or -l^z/L^ 0 (20) 

mm 

$ = (l - a z/L) 1 ^ 2 (21) 

ri ri 



A review of the constants a and a TT has recently been given by 

m H 

2 8 

Busch et al (1973) and fairly good agreement appears to exist 

I 2 9 

for the value of a . Businger et al (1971) suggests a = 15 , 
m m 



o 0 31 

while Paulson (1970), Badgley et al (1972), and Dyer and Kicks 

32 33 

(1970) find a = cr = 16 . Miyake et al (1970) use a = cr =16, 
m K m H 

also used by Pond et al, ^ Phelps and Pond, ^ and Leavitt, 1 1 to obtain 

results from BOMEX data. For neutral conditions (<3? = <J> =1 

\ m H 

29 

Businger et al (1971) find values of = 1. 35 and K = 0. 35 as 

compared to the often used values of Q! = 1, 0 and K = 0.4 , and 

y 

suggest a = 9. 



2. 4 Budget Equations 

If horizontal homogeneity and stationarity are assumed and 

\ 

the effects of humidity on buoyancy are included, the budget 

, . 2 ~2 ~2 ~~Z 

equation for kinetic energy per unit mass, e = u + w + v may 

19 

be written as (Lumley and Panofsky, 1964) 



— b U 
" uw 6~ 



+ [w 0 +0.61 T win ] 
T 

v 



_1 

2 




lb_ 

P 6 z 



pw - c = 0 



( 22 ) 



The non-dimensional form of the budget equation is obtained by 

3 

multiplying Eq. (22) by K z /u and employing Eqs. (15) and (18) so 

that Eq. (22) can be written as 

$ - z/L - D - D - $ = 0 . (23) 

m e p € 



The terms represent the turbulent shear production, buoyant 
production, the flux divergence of turbulent kinetic energy, the flux 
divergence of pressure transport and viscous dissipation 



respectively. 



Three approaches to simplification of Eqs. (22) and (23) have 
been suggested. The first approach, given by Eq. (13), is the 
assumption that production and dissipation of mechanical energy are 

Busch and 
by the sum of 



equal for near neutral conditions. For |z/L | <0.5 
34 



Panofsky (1968) " suggest that dissipation is balancec 

I 

buoyant and mechanical production, and that the flux divergence terms 
may be neglected so that Eq. (23) simplifies to 



<E> - z/E = <£> 

m c £ 



(24) 



35 

McBean, et al (1971) conclude from measurements over land that for 
near neutral conditions (jz/L | < 0. 2) Eq. (24) is a good approxima- 
tion. However, for more unstable conditions, - z/L> 0. 3 , 
dissipation appears to exceed production, and for z/L~ -0. 5 the sum 

§ - z/L - is about 30% of shear production and more than half 

m £ 

33 8 

of buoyant production. Miyake et al (1970) and Pond et al (1971)' 

apply Eq. (24) to calculate momentum flux, however, they do not include 

stability effects on the mean velocity profile as given by Eq. (18). Hicks 
3 6 

and Dyer (1972) include buoyancy effects from Eqs. (18) and (24) to 
calculate the momentum flux by the inertial dissipation technique. 

A recent approach to simplification of the turbulent energy budget 

/ 37 

has been suggested by Wyngaard and Cote (1971) based upon 



\ 

measurements in the surface layer over land. From measurements of 



2 • 2 
we at three different heights the turbulent flux divergence of e was 

estimated and found to approximately balance buoyant production in 

the stability range - 1. 0 £ z /L ^ 0 . Eq. (23) reduces to 



uw 



5 U 

'Oz 



+ e + I 



(25) 



where the imbalance I , is attributed to pressure transport 
divergence. Under stable conditions the imbalance may not be 
substantial and under strongly unstable conditions (z/L = - 2) 
dissipation and total production did balance. Empirically, the 
expression for the imbalance term is written as 



Kz 



u , 



(1-15 z/Lf 1/4 - (1 + 0. 5 | z/L | 2/3 ) 



3/2 



(26) 



where the first term on the right-hand side is the Businger formula 
for dimensionless shear production (Eq, (20), or = 15 ), and the 



m 



second term is the empirical expression determined by Wyngaard and 
_ / 37 

Cote for unstable dissipation data. Using this approach, 
momentum flux corrected for stability effects may be calculated using 
the empirical expression for unstable dissipation written as 



= (1 + 0. 5 | z/L \ 2/i ) 
u_,_ 



3/2 



3 



( 27 ) 



This expression for the dissipation was employed by Stegen et al 
14 

(1973) to calculate momentum flux by the direct dissipation 
technique. 

In recent work employing data from Operation BOMEX, Leavitt 
(1973)^ found that total production approximately balanced dissipation 



in the stability range - 0. 2 £ z/L ^ - 1. 5 , and that 



dimensionless 



turbulent transport as estimated from measurements at two heights was 
only about half of the dimensionless buoyant production as measured by 
Wyngaard and Cote 7 (1971). 

g 

In general the budget of temperature variance (0 ) may be 

27 

written (Monin and Yaglom, 1971) 



be 2 - be 2 — 60 . b 

+ u. : + 2 u-9 t + 



b t 



i b x. 



1 bx. 



£r. ( u i 



9 2 - D. 66 



9 b x. 



- 2D e 'E: 1 + 26R ' + 1 11 + s 9 



(28) 



where repeated indices are summed and u. and x. may be replaced 

i i 

by (u , v , w) and (x , y , z) respectively. Here 0 is the mean 
potential temperature. The temperature variance budget can be 
simplified by invoking the assumptions of stationarity and horizontal 
homogeneity; Eq. (2 8) can then be written as 



w 0 



6 0 
b z 



+ 



1 

2 b z 



w 9 




x e - 29R ' 



2 9c 
C 

P 



- S 



9 



0 (29) 



is defined as 



where the dissipation of temperature variance \ 



0 



x e = 2 d o<Ht> 
1 



2 



( 30 ) 



In a similar manner, the budget of humidity variance 

3 8 

written as (Coantic and Leducq, 1969) 



(q ) may be 
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(31) 



which simplifies under the same assumptions of stationarity and 
horizontal homogeneity to 



wq 



&q 

5 2; 



+ 



1 & 

2 b z 



2 

wq 







(32) 



and the dissipation of humidity variance is defined as 

X' = 2 D ) . (33 ) 

q q 0 x. 

The first three terms of the simplified budget equations for tempera- 
ture and humidity variance (29) and (32), represent production of 
scalar variance, transport of scalar flux divergence, and scalar 
variance dissipation due to the smoothing action of molecular 



diffusivity respectively. The temperature variance budget contains 
additional source or sink terms not included in the humidity budget. 



The term 2 0R' represents the effccvs of heat transfer due to radiation 



on the temperature variance, 



2 0 e 



the effects of viscous dissipation on 



the temperature field due to the internal friction of thfe fluid, and S 



and S represent generalized source or sink terms t 



0 



o account fox' any 

phase changes of water vapor or interactions of the teimperatur e and 
humidity fields due to radiation. 

2 7 39 

Monin axxd Yaglom (1971) emd Plate (1971) suggest *by 



hypothesis that the viscous internal heating term 20C/C is insignificant 

P 

40 

and may be neglected. Friehe (1973) suggests that this term may not 

be negligible and that there may be a positive correlation of 0 C due to 

an increase in 0 produced by local, rapid adiabatic heating by viscous 

dissipation. Although a mean valxxe for viscous dissipation £ may 

^~2 

be estimated from the mean square of the velocity derivative (-p ) 

dT^ 

(see Eq. (13)), a positive correlation of 0 (pp) may not accurately 

reflect the actual correlation between fluctuating temperature and 



the instantaneous viscous dissipation Qf ; (— ) is only one component 
of dissipation, used to calculate £ by the assumption of isotropy. 



In general, the effect of radiative heat transfer is to reduce the 

41 

value of the mean square temiperature fluctuation (Goody (1956), 

42 39 

Townsend (1958), Plate (1971) ). For internal radiative transfer 



occurs at 



42 2 

Townsend indicates that maximum destruction of 0 



very small scales of the temperature fluctuations since the most 

pronounced radiative transfer occurs over the shortest path interval. 

For atmospheric conditions with high humidity content and strong 

JO 



radiative transfer, Phelps and Pond (1971) suggest 
absorption of long wave radiation occurs at larger sc 



that more 
lies resulting in 



a more homogeneous large-scale temperature field and suppression of 

39 

low frequency temperature fluctuations. Plate also indicates that 
partial absorption of long-wave radiation by atmospheric water vapor 
could cause a change of state**, (conver sion to sensible heat) which 
would influence the temperature field. It is not clear whether this 
absorption and change of state of water vapor would act as a source or 
sink in the temperature and humidity variance budgets, and could be 
included in the generalized source terms. 



Coantic and Seguin (1971) have shown that radiative flux 
divergence due to absorption of inf a red radiation (long-wave) by water 
vapor and carbon dioxide always acts in such a way to increase the 
absolute value of the turbulent heat flux. They conclude that the constant 
heat flux hypothesis for the surface layer may not be valid, especially 
for conditions of moderate wind and high moisture content, even for 
heights above 1 meter. 



The generalized source term may also include the effects of ocean 



